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This is the first of series of papers in which we investigate stability of the spherically symmetric 
space-time with de Sitter center. Geometry, asymptotically Schwarzschild for large r and asymptot- 
ically de Sitter as r — > 0, describes a vacuum nonsingular black hole for m > m^r and particle-like 
self-gravitating structure for m < mcr where a critical value rricr depends on the scale of the sym- 
metry restoration to de Sitter group in the origin. In this paper we address the question of stability 
of a vacuum non-singular black hole with de Sitter center to external perturbations. We specify first 
two types of geometries with and without changes of topology. Then we derive the general equations 
for an arbitrary density profile and show that in the whole range of the mass parameter m objects 
described by geometries with de Sitter center remain stable under axial perturbations. In the case 
of the polar perturbations we find criteria of stability and study in detail the case of the density 
profile p(r) = poe~ 



where po is the density of de Sitter vacuum at the center, ro = 



is de Sitter radius and rg is the Schwarzschild radius. 



PACS numbers; 04.70.Bw, 04.20.Dw 
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I. INTRODUCTION 

The idea of replacing of a Schwarzschild singularity 
with de Sitter vacuum goes back to 1965 papers of 
Sakharov [1] who considered p — —p as the equation of 
state for superhigh density and of Gliner who interpreted 
p ~ as corresponding to a vacuum and suggested that 
it could be a final state in a gravitational collapse [2] . 

In 1968 Bardeen presented the spherically symmetric 
metric of the same form as Schwarzschild and Reissner- 
Nordstrom metric, describing a non-singular black hole 
(BH) without specifying the behavior at the center [3]. 
The very important point was noted in [3] for the first 
time: that the considered space-time exhibits the smooth 
changes of topology. 

Direct matching of Schwarzschild metric to de Sitter 
metric within a short transitional space-like layer of the 
Planckian depth [4-8] results in metrics typically with a 
jump at the junction surface. 

The situation with transition to de Sitter as r — s- 0, 
has been analyzed in 1988 by Poisson and Israel who 
found necessary to introduce a transitional layer of " non- 
inflationary material" of uncertain depth at the charac- 
teristic scale {r^TgY/^ (ro is de Sitter radius, and rg is 
the Schwarzschild radius), where geometry can be self- 
regulatory and describable semiclassically down a few 
Planckian radii by the Einstein equations with a source 
term representing vacuum polarization effects [9]. 

Generic properties of " noninflationary material" have 
been considered in 1990 in Ref. [10]. For a smooth de 
Sitter-Schwarzschild transition a source term satisfies [10] 



and the equation of state, following from = 0, is 



-p; Pi- 



(1.2) 



Here k = SttG (we adopted c = 1 for simplicity), 
p{r) = Tl is the energy density, Pr{r) = — is the ra- 
dial pressure, and p^_{r) = —Tg = —T^ is the tangential 
pressure for anisotropic perfect fiuid [11]. 

The stress-energy tensor with the algebraic structure 
(1.1) has an infinite set of comoving reference frames and 
is identified therefore as describing a spherically symmet- 
ric vacuum [10], invariant under boosts in the radial di- 
rection and defined by the symmetry of its stress-energy 
tensor (for review [12-16]). 

The exact analytical solution was formd in 1990 for the 
case of the density profile [10] 



p(r) =poe rl^S/npo; 



rg — 2Gm 



(1.3) 



t:- n 



(1.1) 



which describes a vacuum asymptotically de Sitter as r — s- 
in a simple semiclassical model for vacuum polarization 
in the spherically symmetric gravitational field [17]. 

In 1991 Morgan has considered a black hole in a simple 
model for quantum gravity with quantum effects repre- 
sented by an upper cutoff on the curvature, and obtained 
de Sitter-like past and future cores replacing singularities 
[18]. In 1992 Strominger demonstrated the possibility of 
natural, not ad hoc, arising of de Sitter core inside a black 
hole in the model of two-dimensional dilaton gravity con- 
formally coupled to N scalar fields [19]. 

In 1996 it was shown that in the course of Hawking 
evaporation a vacuum nonsingular black hole evolves to- 
wards a self-gravitating particle-like vacuum structure 
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without horizons [17], kind of gravitational vacuum soli- 
ton called G-lump [20]. The form of temperature- mass 
diagram is generic for de Sitter-Schwarzschild black hole 
[17] and dictated by the Schwarzschild asymptotic and 
by the existence of two horizons - when decreasing dur- 
ing evaporation mass reaches a certain critical value rricr 
evaporation stops [17,21]. 

In 1997 in Ref. [22] it was shown that in a large class 
of space-times that satisfy the Weak Energy Condition 
(WEC), the existence of a regular black hole requires 
topology change. Bardeen metric and the metric gener- 
ated by the density profile (1.3) belong to this class. 

In 2000 studying the quantum gravitational effects by 
the effective average action with the running Newton 
constant, and improving Schwarzschild black hole with 
renormalization group, Bonnano and Renter [23] have 
constructed nonsingular black hole metric and confirmed 
the results of [17] concerning the form of temperature- 
mass diagram and the fundamental fact that evaporation 
stops when the mass approaches the critical value rricr- 

Also in 2000 the regular BH solution with a charged de 
Sitter core has been considered by Kao [24] with using the 
density profile (1.3) for distribution of a charged material. 

In the same year 2000 regular magnetic black hole and 
monopole solutions arc found by Bronnikov [25] in Non- 
linear Electrodynamics (NED) coupled to gravity with 
the stress-energy tensor of the structure (1.1). 

Existence of regular electrically charged structures in 
nonlinear electrodynamics coupled to general relativity 
was proved recently in Ref. [26], where it was shown 
that in NED coupled to GR and satisfying WEC, reg- 
ular charged structures must have de Sitter center. 

In 2001 the non-singular quasi-black-hole model rep- 
resenting a compact object without horizons, was con- 
structed by Mazur and Mottola [27] by extending the con- 
cept of Bose-Einstein condensation to gravitational sys- 
tems. An interior de Sitter condensate phase is matched 
to an exterior Schwarzschild geometry of arbitrary mass 
through a phase boundary of a small but finite thickness 
with equation of state p = p. 

In 2002 nonsingular BH solution was found by Nashed 
[28] as a general solution of Moller tetrad theory of gravi- 
tation by assuming the same specific form of the vacuum 
stress-energy tensor as in Ref. [10] with the density pro- 
file (1.3). Later it was extended to the case of teleparallcl 
theory of gravitation [29] . Stability condition of geodesic 
motion in the field of vacuum nonsingular black hole de- 
scribed by the regular analytic solution [10] with the den- 
sity profile (1.3), has been considered in [30]. 

Model-independent analysis of the Einstein spheri- 
cally symmetric minimally coupled equations has shown 
[20,31] which geometry they can describe if certain gen- 
eral requirements are satisfied: (a) regularity of density; 
(b) finiteness of the ADM mass; (c) Dominant Energy 
Condition (DEC) for Tj^i/. These conditions lead to exis- 
tence of regular structures with de Sitter center including 



regular black holes without topological changes. The ex- 
ample of such a case is the exact analytic solution [26] 
describing in certain mass range a regular charged black 
hole with de Sitter center. 

The condition (c) can be loosed to (c2): weak energy 
condition for T^^ and regularity of pressures [31]. WEC 
which is contained in DEC, in both cases is needed for 
de Sitter asymptotic at the center. 

The requirements (a)-(c) either (a)-(c2) define the fam- 
ily of asymptotically flat solutions with the regular center 
which includes the class of metrics asymptotically de Sit- 
ter as r — > 0. A source term connects de Sitter vacuum in 
the origin with the Minkowskli vacuum at infinity. Space- 
time symmetry changes smoothly from de Sitter group 
at the center to the Poincare group at infinity through 
the radial boosts in between, and the standard formula 
for the ADM mass relates it to both de Sitter vacuum 
trapped inside an object and smooth breaking of space- 
time symmetry [20]. 

Cases (c)-(c2) differ by behavior of the curvature scalar 
R. In the case (c) it is non-negative which evidences 
the existence of regular black holes without topological 
changes. So, the class of metrics with de Sitter cen- 
ter includes two subclasses with and without topological 
changes. 

In this paper we specify conditions of existence of two 
types of geometries with the de Sitter center and investi- 
gate stability of configurations described by these geome- 
tries, by studying perturbations in geometry via Einstein 
equations linearized about the unperturbed space-time. 
Results are valid for geometries of both types. 

This paper is organized as follows. In Sect. 2 we out- 
line the conditions of existence and basic properties of 
spherically symmetric geometries with de Sitter center. 
In Sect. 3 we introduce the basic equations describing ax- 
ially symmetric time-dependent perturbations of a spher- 
ically symmetric system with dc Sitter center. In Sect. 4 
we prove stability of such a system to axial perturbations. 
In Sect. 5 we analyze the case of polar perturbations and 
derive criteria of stability to these perturbations. In Sect. 
6 we apply the results to the case of the vacuum nonsin- 
gular black hole with the density profile (1.3). Section 7 
contains summary and discussion. 

II. SPHERICALLY SYMMETRIC SPACE-TIME 
WITH DE SITTER CENTER 

A static spherically symmetric line element can be 
written in the form [11] 

ds'^ = e'^^^'Ut^ - e'^^^Ur'^ - r^dft'^ (2.1) 

where dfl'^ is the metric of a unit 2-sphere. Integration 
of the Einstein equations gives 

e—(r) ^ 1 _ 2GM(r) _ M{r) = 4n £ p{x)x'^dx (2.2) 
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whose asymptotic for large r is e = 1 — 2Gm/r, with 
the mass parameter 



m 



= 47r / 
Jo 



p{r)r^dr 



(2.3) 



Requirement of regularity of density, po = p{r — > 0) < oo, 
loads to behavior of mass function M(r) ^ as r — > 
and thus i/(0) = 0. 

To outline the conditions of existence of spherically 
symmetric space-time with de Sitter center, we need the 
Oppenheimer equation [32] 



T;=pr+p- 



1 e" 



(2.4) 



The dominant energy condition T^" > |T"''| for each 
a,b = 1,2,3, holds if and only if [33] p > 0; p + Pk ^ 
0; P — Pfc > 0; k = 1,2,3. It includes the weak energy 
condition which implies p > 0; p+pk > 0. Together with 
the condition of regularity of density, DEC (via pk < p) 
leads to ju' + I/' = as r ^ [20]. 

The same result can be achieved by requirement of 
regularity of pressure (the subclass satisfying (a-c2)). 

In the limit r ^ oo the condition of finiteness of the 
mass (2.3) requires density profile p{r) to vanish at in- 
finity quicker than r~^. In the case (c) the dominant 
energy condition requires pressures to vanish as r ^ oo. 
Then p' = and p =const at infinity. Rescaling the 
time coordinate allows one to put the standard bound- 
ary condition — » as r ^ oo which ensures asymptotic 
flatness needed to identify (2.3) as the ADM mass [34]. 

The same result can be achieved in the case (c2) by 
postulating regularity of pressures including vanishing of 
Pr at infinity sufficient to get p' = needed for asymp- 
totic flatness. 

The weak energy condition requires p' + u' > 0. The 
function p + u is growing from = ^(0) at r = to /x = 
at r ^ DC', which gives /i(0) < [20]. 

The range of family parameter ^(0) includes p{0) = 0. 
In this case the function j/(r) -I- p{r) is zero at r = and 
at r oo, its derivative is non- negative (by WEC via 
P+Pk > 0), it follows that i/(r) = —p{r) everywhere. 

A source term for this class of metrics corresponds to 
anisotropic perfect fluid which satisfies the r— dependent 
equation of state (1.2), and the weak energy condition 
p± + p > demands monotonic decreasing of a density 
profile, p' <0 [20]. 

Behavior at r -> is dictated by the WEC [20]. The 
equation of state near the center becomes p = —p , which 
gives de Sitter asymptotic as r — > 



dr^ 



\ 



A = Kpo 



(2.5) 



(2.6) 



where pQ = p{r 0) and A is the cosmological constant 
which appeared at the origin although was not present 
in the basic equations. 

Requirements (a-c) either (a-c2) lead thus to the exis- 
tence of the class of metrics 



dr^ 



ds'^ = g{r)dt^ - ^ - r'^dn^ 



9{r) 



9{r) 



Rg{r) = 2GM{r) 



M{r 



) = 47r / p{x)a 
Jo 



'dx 



(2.7) 



(2.8) 



(2.9) 



which are asymptotically de Sitter as r — > 0, and asymp- 
totically Schwarzschild at large r 



ds' 



1-^ 

r 



dr^ 



1- ^ 



-r^dn^; rg = 2Gm (2.10) 



The weak energy condition defines the form of the metric 
function g{r). In the region r > it has only minimum 
and the geometry can have not more than two horizons: 
a black hole horizon r_|_ and an internal horizon R_ [20]. 

The scalar curvature R, proportional to the trace of 
stress-energy tensor T, is proportional to p — p± for ge- 
ometries satisfying (1.2), so that conditions (a-c) and (a- 
c2) distinguish two types of geometries. In the case (a- 
c) satisfying DEC requirement, scalar curvature remains 
non-negative, since DEC requires p — Pk > 0. The sub- 
class satisfying (a-c) does not exhibit changes of topology 
by virtue of DEC and can be specified as DEC-subclass. 
Dominant energy condition requires that each principal 
pressure does not exceed the density which guarantees 
that speed of sound can not exceed speed of light. In 
nonlinear electrodynamics coupled to gravity, photons 
do not follow null geodesies of background geometry but 
propagate along null geodesies of an effective geometry 
[35], and propagation of photons resembles propagation 
inside a dielectric medium [36] . In the case of the regular 
NED structure satisfying DEC [26] , it allows one to avoid 
problems with speed of sound exceeding speed of light. 

In the case (a-c2) scalar curvature R{r) changes 
sign somewhere and geometry experiences topological 
changes. This subclass satisfying only weak energy con- 
dition (needed in both cases for de Sitter behavior at the 
center) can be specified as WEC-subclass. 

The case of the density profile (1.3) belongs to WEC- 
subclass satisfying (c2). The metric function and the 
mass function are given by [10] 



7(r) 



1-^ 
r 



3/2 



M(r) 



1- 



-r^/r 



(2.11) 
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Dominant Energy Condition is not satisfied so that a 
surface of zero scalar curvature exists at which R{r) = 0. 
Zero curvature surface r = Tg is shown in fig.l together 
with two horizons (a black hole event horizon r+ and 
an internal horizon r_ ) , and the characteristic surface of 
any geometry with de Sitter center: a zero-gravity sur- 
face r = Vc beyond which the strong energy condition of 
singularities theorems [33], is violated (zero-gravity sur- 
face is defined by 2p + rp' = [17]). 




FIG. 1. Chaxacteristic surfaces of a spherically symmetric 
space-time of WEC type with de Sitter center. 

Two horizons come together at the value of a mass 
parameter rricr, which puts a lower limit on a black hole 
mass (sec fig. 2). For the case of a density profile (1.3) 
the critical mass is given by [17] 



(2.12) 



1 .5 
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FIG. 2. Metric function g{r) for de Sitter-Schwarzschild 
configurations. Mass m is normalized to rricr- 

Temperature-mass diagram is shown in fig. 3. Its form 
does not depend on particular choice of a density pro- 
file. Temperature drops to zero at m = rricr, while the 
Schwarzschild asymptotic requires T_|_ ^ as m ^ oo. 
As a result the temperature-mass diagram should have 
a maximum between mcr and m, ^ oo [17]. In a maxi- 
mum, at m = mcr2, a specific heat is broken and changes 
sign testifying to a second-order phase transition in the 
course of Hawking evaporation [21]. 




FIG. 3. Temperature-mass diagram for a vacuum nonsin- 
gular black hole with de Sitter center. 

For m > rricr de Sitter-Schwarzschild geometry de- 
scribes the vacuum nonsingular black hole, and global 
structure of space-time shown in fig. 4 [17], contains an 
infinite sequence of black and white holes whose future 
and past singularities arc replaced with regular cores TZC 
asymptotically de Sitter as r ^ [17]. 




FIG. 4. Global structure of space-time of a vacuum non- 
singular black hole with de Sitter center [17]. 
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III. BASIC EQUATIONS FOR PERTURBATIONS 

The perturbations of a spherically symmetric system 
arc on essence time-dependent axially symmetric modes; 
the reason is the absence of a preferred axes in a spheri- 
cally symmetric background [37]. 

They are described by the line element [37] 

ds^ = e^^dt^ - e^^(# - ojdt - q2dr - q^dOf 

-e'^'^'idrf -e'^'^^idOf, (3.1) 

in which metric functions i^, i/;, /i2, w, g2, 93 are func- 
tions only of t, r, 6. They satisfy the Einstein equations 

Gij = Rij — —gijR = ~KTij, (3-2) 

where the Ricci tensor is defined by Rij = g''''Rikji- 
Non-zero components of stress-energy tensor read 

Trr = e'^i'^pr, Tee = e'^'^^pe- (3.3) 

where Pr,Pe,Pci> are the principal pressures. 

We obtain the relevant perturbation equations by lin- 
earizing the field equations around the spherically sym- 
metric solution with de Sitter center. This solution con- 
sidered as a special case of the line element (3.1) with 

M2 = -fir); V' = ln(rsin(6')); ^3 = ln(r), (3.4) 

has the form (2.7) with 




The particular solution (3.5) is specified by the choice of 
the constant Ci which we choose in such a way to have 
unperturbed metric given by (2.8)-(2.9). 

Our task is to investigate stability of the spherically 
symmetric system with de Sitter center to external per- 
turbations in general case of a regular density profile p{r). 

The class of metrics with de Sitter center and a source 
term of structure (1.1), is extended to the case of non- 
zero cosmological constant (A < A) at infinity [38] corre- 
sponding to extension of the Einstein cosmological term 
Agik to an r— dependent second rank symmetric tensor 

A,fe = KTjfe (3.6) 

with the algebraic structure (1.1), connecting smoothly 
two de Sitter vacua with different values of a cosmolog- 
ical constant [39]. In this approach a constant scalar A 
associated with a vacuum density A = Kpvac, becomes a 
tensor component A* associated explicitly with a density 



component of a perfect fiuid tensor whose vacuum prop- 
erties follow from its symmetry (1.1) and whose variabil- 
ity follows from the Bianchi identities [39,20]. 

Here we investigate stability for the particular case 
when A = and spherically-symmetric space-time with 
de Sitter center is asymptotically flat. 

Since an anisotropic fluid with the stress-energy ten- 
sor of type (1.1) admits identifying it as a vacuum-like 
medium associated with a time-evolving and spatially 
inhomogeneous cosmological term [10,39,20,40], we can 
write the Einstein equations in the form 

Gik+Aik = (3.7) 

(for discussion of where to put cosmological term see 
[41,20]). Then the quantities p,Pk are treated as corre- 
sponding (in one-to-one way) components of the variable 
cosmological term A* = np,h\ = —upk [39]. 

Since we apply the approach of studying direct pertur- 
bations of geometry via Einstein equations, we consider 
behavior of small perturbations for both the metric ten- 
sor and a stress-energy tensor associated with A^^. 

A general perturbation of a background geometry will 
result in ui,q2,q3 becoming small quantities of the flrst 

order and the functions v, p2, PStiP and p, pk experiencing 
small increments 5v, 5p2, Spa, dtp and Sp, Spk- 

The perturbations leading to non-vanishing values of 
LU,q2 and qs induce a dragging of the inertial frame and 
impart a rotation, for this reason they are called axial 
perturbations [37]. 

Perturbations which do impart no rotation are called 
polar perturbations [37]. In the considered case they lead 
to increments in i',p2,P3,i' and p,Pk- 

The equations governing the axial and the polar per- 
turbations decouple. 

Axial perturbations are governed by equations [37] 

Rr4> = Re4> = (3-8) 
The equations governing the polar perturbations read 

-Rtr = + P3),rt + i^,r{i^ ^ P2) ,t 
+M3,r(M3 - P2),t - i^,r(V' + l^s) ,t = 0, (3.9) 

-Rte = {tp + P2),et + i',e{tp - P3),t 
+P2,e{P2 - P3),t - i^,e{'^ + P2),t = 0, (3.10) 

-Rre = + v)^re + V',r(V' " P2) ,B 
+V^r{v - P2),e - P3,r{-4' + f) ,0 = 0, (3.11) 
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Gtt=e ^''^[(V'+M3),r-r+V',r(V'-M2+M3),r+M3,r(M3-M2),r] Nowwe take 



+e [{iIj + ij.2),ee + '^,e{4' + M2 - /U3),9 + 1^2,9(1^2 - 1^3), e] 



-e-'"[V',t(M2 + /X3),t + M2,tM3,t] = -Ke^V, (3.12) 

G^,/, = e~^''^[(!/+/i3),rr+!^,r(z^-M2+M3),r+/"3,r(M3-/"2),r] 



-e ^"[(^2 + M3),<t + A'2,t(M2 + M3 - f),t + l^3,t{l^3 " ^) ,t\ 



= (3.13) 



Grr = e ^'^^ [V',r(i^ + M3),r + i^,rA*3,r] 



-e ^""[(V) + /U3),« + V',t(V' + M3 - z^),t + M3,t(M3 - z^),t] 



(3.14) 



Gee = e ^^''[(■i/' + J^),rr + V',r('i/' + I^-A*2),r + I^,r(l'-A<2),r] 



+V,t(V' + M2 - '^),t + /"2,t(/U2 - !^),t] = Ke'^'^'pe- (3.15) 

We perturb equations (3.8)-(3.15) up to the first order, 
and in equations (3.12)-(3.15) we disturb both left and 
right hand sides. As a result we obtain the linear system 
of 7 partial differential equations for the polar perturba- 
tions, and the linear system of 2 equations for the axial 
perturbations. 



IV. AXIAL PERTURBATIONS 

Axial perturbations corresponds to appearing of 
nonzero values w, (72, 93 which vanish for unperturbed sys- 
tem. 

They are governed by the Einstein equations (3.8). 
This gives 2 equations for 3 functions which read 

eMr) 

[sin e{q2,e - Q3,r)],e = -(w,r - q2,t),t (4.1) 



o . 3 

r'' sin 



[r^e^-'^^\q2,e - qs,r)],r = {uJ,e - q3,t),t (4.2) 



and similarly for (72,93; in what follows we retain the 
same symbols for the amplitudes of the perturbations 
which satisfy equations 



sin'^ t 



sin^ 0{q2,e - 93,r 



r'e2^W(92,e-g3,r) 



= -{iauj,r + (T^q2) (4.3) 



iato^e + a^qs (4.4) 



Expressing ui^r from (4.3) and ui^g from (4.4), differenti- 
ating and equating w^rO = ^,er we get one equation 



(4.5) 



dr V dr 
Here 

Q(r, 9) = e^^r' sin^ e[ci2.B - gg.r) (4.6) 



dr \ dr ) 



d ( 1 de\ A ^ „ 



For the Schwarzschild metric, eq.(4.5) coincides with the 
analogous Chandrasekhar equation ( [37], Ch.4, eq.(18)). 
Separating variables by Q{r,9) = R{r)Q{6), we get 

.V4(^?)-Aiii + .^« = (4.7) 

M\^xFellB ) ' sin^"~" ^^'^^ 
Solutions to (4.8) are Gegenbauer polynomials 

Qm = G;";|(^) = {Pifie - PifiCtgO) sin^ 6 (4.9) 
which gives 

A; = (Z + 2)(Z-1); Z = 2,3,... (4.10) 

General solution can be written in the form 

00 

Q{r,e) = Y,Ri{r)ei{e) (4.11) 



1=2 



Equation for Ri (r) reads 



dr V dr 



-^{l + 2){l-l)Ri + afRi = 

(4.12) 

In "tortoise" coordinate = J dr/g{r), we get 
Schrodinger equation for Zi{r^) = r~^Ri{rt,) 



(4.13) 
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where the potential is given by 

Vi{r) = ^ {li^ + 26^" - 2rure^'') (4.14) 

with 

fi^ = {l + 2){l-l) (4.15) 

For the Schwarzschild geometry (4.14) coincides with 
the Rcggc-Whcclcr potential ( [37], Ch.4, eq.(28)). 
With the asymptotic behavior of (4.14) 

Vi — — as r ^ r* ^ oo 

Vi const e9'(''+)''* as r* ^ -oo (4.16) 

solutions of (4.13) have asymptotic e^"^""* as r* ±oo 
as in the case of the Schwarzschild geometry [37]. For 
real a they describe propagation of ingoing and outgoing 
waves through one-dimensional potential barrier, so that 
we have to look for solutions to (4.13), which satisfy the 
boundary conditions [37] 

Zi ^ e+^""^* + i?,(az)e-^""^* (n ^ +oo) 

Zi ^ Ti{ai)e+''''^' (n ^ -oo) (4.17) 

These boundary conditions tell us that each I component 
of an incident wave of the unit amplitude coming from 
infinity gives rise to a reflected wave of amplitude Ri{<Ji) 
at infinity and a transmitted wave of amplitude Ti{ai) at 
the horizon [37]. 

If such solutions exist only for real values of the time 
parameter a; and form complete basic set, then any 
smooth initial perturbation defined at the finite interval 
of (with compact support), can be expanded on these 
functions, and since dependence of perturbations from 
time coordinate has the form exp(z(Tji), this is followed 
by stability of geometry in question. 

In terms of the metric function g{r) 

Vi{r) = ^-^{„^ + 2g-rg') (4.18) 
It is easily seen that Vi (r) is constrained by the function 

Vi{r) > ^[l + g+\'^pr''y 

which is certainly positive and presents the value of the 
potential for the mode 1 = 2. 

We see that axial perturbations are governed by one- 
dimensional wave equation (4.13) with a non-negative po- 
tential. In terms of a one-dimensional Schrodinger equa- 
tion, an unstable mode exists if a potential has a bound 



state, which corresponds to the negative eigenvalue af . 
In the case of non-negative potential the system obeys the 
theorem [42] which guarantees the absence of negative 
eigenvalues in the standard one-dimensional Schrodinger 
equation with the non-negative potential. The absence of 
negative eigenvalues in the spectrum of (4.13) guarantees 
the absence of exponentially growing modes. 

For the case when the potential is real, smooth and 
short-range, standard theorems of quantum mechanics 
guarantee that cigcnfimctions of (4.13) form a complete 
set and any square-integrable state function can be ex- 
panded on them [37]. 

We have however to be careful about the behavior of 
the solution to (4.13) in the extremal regime near the 
double horizon r± which satisfies g{r±) = 0; g'{r±) = 0. 

To study the extreme case we introduce dimensionless 
variable x by normalizing the variable r to the character- 
istic scale {r^TgY^^. 

Then the equation for the function Ri reads 

x^g^R'{ + 2x(^-gyR:^ + {a^x^-^i'g)Rl = (4.19) 

where oi = ailr'^rgY^^. In what follows we retain the 
notation ai keeping in mind that it is multiplied by 
(rgr^)^/^. (A multiplier is not essential since in studying 
stability we are interested only in the sign of af.) 

Near the double horizon x±, the metric function is 
g{x) = g"{x±){x — x±f 12 + .... Introducing the vari- 
able z = X — x± and the notation 7 = g"{x±)/2, we get 
the limiting equation 

2 

z^Ri^,, + 2z''Ri,, + ^Ri = (4.20) 

General solution to (4.20) (found by taking Ri = 
ivi/z), reads [44] 

Riiz) = Cucos(^'^^ + C2isin(^'^^ (4.21) 

(with = o-f /j"^), and tells us that axial perturbations 
are restricted near double horizon. 

In case without horizons asymptotic behavior of the 
potential at infinity and near zero is Vi = l{l + l)/r'^. 
This is behavior of centrifugal part of radial Schrodinger 
operator in the spherically symmetric case. 

We see that potential for axial perturbations is smooth, 
short-range and positive for all types of configurations de- 
scribed by spherically symmetric geometry with de Sitter 
center. We conclude that geometry with de Sitter center 
is stable to axial perturbations. 
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V. POLAR PERTURBATIONS 

A. General equations 

Linearizing the system (3.9)-(3.15) about the back- 
ground geometry, we get the system describing polar per- 
turbations 



0, (5.1) 



m + Sf^2),e + ctg{9){Sij - (5/X3)],t = 0, (5.2) 



u_r - - iv^r + - )<5/X2,e = 0, (5.3) 



— e 



2v 



{Sip + 5/U3),rr + 



--^/X2,r--2( -I^,r- + -2 



r 



„2l/ 



+26^3 + 5^2,66 + ctg{0)6ii2fi] = K.Sp, (5.4) 



{61/ + 5fJb3)^rT + ( 3z^,r- + ^ V^^.r- + 2^ I^,r + ^ ) <5A*3,r 



l',r + - j ^M2,r — 2 iv^rr + 2v'^^ + -P^r ) '^M2 



+^{5i^ + 5fi2),ee - e '^''{Sij.2 + Sfi3)^u = ^Sp^,, (5.5) 



"2^ 












_r 







„2!^ 



+^[6^,ee+ctg{0){2Stp-Sii3),0+2Sn3+Si',ee+ctg{0)Siy,e] 



-e "^"{Sip + diJ,3)^tt = nSpr, (5.6) 

((5z^ -h Slp)^rr + (^'^,r + -^Sv^r + 2^,r + '^V'.r 



U^r + ^jSH2,r - 2 (l',rr + 2u'^^ + ^V^r ] 5^2 



+^ctg{e){Su + Sfj,2),e-e~'^''{Sfj,2 + S^),tt = k5p0. (5.7) 

The system (5.1)-(5.7) is the system of 7 Hnear par- 
tial differential equations of the first order for 8 quanti- 
ties: 4 small perturbations of metric tensor and 4 small 
perturbations of stress-energy tensor (which is in consid- 
ered case can be associated with a variable cosmological 
term) whose unperturbed components are related by the 
equation of state, in our case (1.2). To investigate this 
system we should make an assumption concerning per- 
turbation of Pr valid for the case of small perturbations. 
Since for the background geometry we have Pr = —p, i.e. 
Pr = Pr{p), we can assume (see, e.g., [43]) 



which results in 



dpr 

Spr = -J- Op. 

dp 



Spr = —5 p. 



(5.8) 



(5.9) 



The possibility to connect perturbations 5pr and 5p is 
implied by our system which contains 7 equations for 
8 functions. The relation (5.9) is valid only for small 
perturbations, since only in this case the relation (5.8) 
is valid. So, if we prove that the system is stable, i.e. 
growing perturbation modes are absent, this will justify 
the validity of (5.9). 

Taking into account (5.4) and (5.6), the equation (5.9) 
can be written in the form 



2 2 2 

-{5i> + ()>3),rr {5i' + 5n3)^r + -^V r + -(5/Lt2,T 

r r r 



5v^ee + ctg{9)6ufi — SiJ,2,ee — ctg{9)5ii2,{ 

-e-2'^(J?A + (5/i3),tt = 0. 



(5.10) 



In this way we obtain the system of 7 equations for 7 
unknown functions which splits into uniform system of 
4 linear partial differential equations (5.1), (5.2), (5.3), 
(5.10) for 4 small perturbations of the metric tensor, 
Su{r,6,t),6^,2{r,6,t),5iJ,3{r,6,t),6tp{r,6,t); and 3 linear 
algebraic equations (5.6), (5.5), (5.7), determining ex- 
pressions for Spr, Sp^, Spe through expressions for metric 
perturbations. 

The problem ultimately reduces to investigation of the 
uniform linear system (5.1), (5.2), (5.3), (5.10). 

Following Chandrasekhar [37] we assume the time de- 
pendence e"^* which corresponds to the Fourier analysis 
of perturbations. The variables r and 6 are separated by 
the Friedman substitutions [37]. 

As a result we present perturbations as series 



-|-oo 



(5.11) 



1=2 
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+ 00 / 1\ 1 

Sii2{r,e,t) ^J2Li{r)Pi{cose)e"''\ (5.12) (^iW+^iW),r- (^;^,r--J (^/(r) + L;(r)) + -L;(r) = 0; 

1=2 



+ 00 



1=2 



+Vi{r)PiMcose)]e 



iTit 



(5.13) 



{Ti{r)-Vi{r)),r-[ v,r+^Li{r)+Ni^r{r) + L,r-\ ]Ni{r) = 0; 



„2v 



Xi{r) + Li{r) ] +-[Ni{r) + 2Xi{r) + ZLi{r. 



+ 00 



6ij{r,e,t)^Y.^Ti{r)Pi{cos0) 



+Vi{r)Pifi{cose)ctgey"'\ (5.14) 



+ 00 

5p(r,e,t) =Y,Ci{r)Pi{cos9)e'^'\ (5.15) 



Ni{r)-Li{r)] - e'^^af { Xi{r) + Li{r) ] =0 



1(1 + I) 
2r2 



With using these relations wc transform uhimatcly our 
starting system to the system of 3 differential equations 
in the normal form for the functions Ni{r),Li{r),Xi{r) 

Ni,r = {n + l)alNi + {v^r + hi 



-{n + l)al + afcl)Li + (61 + afcl)Xi, (5.20a) 



Spr{r,e,t) = Y,Di{r)Pi{cose)e'"'\ (5.16) 



1=2 



Li,r = {Kr-^ + {n + l)al\Ni + (-^ + bl 



+ 00 



5p^{r,0,t) =J2[Eiir)Pi{cose) 



-(n + l)al + afclj Li + (61 + afcl))Xi, (5.206) 



+Hi{r)Pi,e0{cos9)y^'*, (5.17) 



+ 00 



Sp0{r,e,t) = ^[Ei{r)Pi{cose) 



+Hi{r)Pi,e{cose)ctge]e''''K (5.18) 

Let us introduce the function Xi which will be useful 
in our further reductions 



Xi{r)=nVi{r), (5.19) 



where 



n=^^^-^-l; Z = 2,3,...; n = 2,5,9,.... 



Using the properties of the Legendrc polynomials 
{sinePifi)fi + l{l + l)sin0Pi{cos0) = 0, 



N, 



+ ^p^r - - 61 + (n + l)al - afcl^ Li 

+ ^ly^^ - 1 - 61 - afcl^ Xi , (5.20c) 



where 



-2v(r) 



al(r) = ; 61 (r) = —rv^rr — i'^r'^ cl(r) = re 

(5.21) 

and 4 equations which define amplitudes Di{r), Ei{r), 
Hi{r), and Ci{r) through solutions of (5.20) 



A(r) = 



\^l,r- + \ ]{Xl+Ll),r 



PlM + PlfiCtgO = -l{l + l)Pi{cos0), 

we get from Eqs (5.1)-(5.3), (5.10) after some algebra, 
the following relations between amplitudes 

Ti{r)=Vi{r)-Li{ry, 



2 1 , . 



:[Xi - nil + {n + l)Ni] 



afe-^-iXi+Li) 



(5.22) 



9 



Ei{r) = 



N1-L1 + -Xi 
n 



1 , 1 (n + 1) 
25 



X xg 



Ni + 



1 g 1 X I g \ X g 



25 



2 5 



+ ( 3i/,^ + J)^'.'- + 2('^,r + J j (^^( - Li 



+cr, e —A; 



(5.23) 



HK 



Xl^rr + 2[ + - )-^i,T 



+-{Ni + Li)+ate-'-Xi 



(5.24) 



CKr) = -A(r), (5.25) 
Let us now introduce the dimcnsionless variables 



X = — ; p ^ — ; where = r^r 
ri po 

and the characteristic parameter 



(5.26) 



a = — 



(5.27) 

In these notations the unperturbed solution (2.7) reads 
ds'^ = g{x)dt^ - 



-2 x^dQ^ 



9{x) 



9{x) 



1 - ^Mi^l. Mix) =3 r piq)q^dq (5.28) 
X Jo 



In terms of g{x) our basic system (5.20) takes the form 

2 



Li + 


[II 











xg 



(n + 1) 2 X 

^ '^l~o 

xg 5 



"ig' 1 I (n+l) 
.2 fif a; a;ff 

a;5f 



2g^A\g) 2g 



xg 



X, 



(5.29a) 



Ni- 



1 a; / 5' \ a; g" 



a; 4 Vffy 2 g 



27 + "'7 



Xi, 
(5.296) 



(n+1) 2 a; 

H fT; ^ 

a;^ g^ 



1 g' 1 X / g'\ X g" ^x 

^'+[27-^-4 +27""'? 

(5.29c) 

This system can be transformed to the equivalent form 
xg'Ni,^ = (n + \)gNi + g (^|g' - (n + 1)^ U 



Xr. 



W[\{g'f-\gg" + o]\X, 



(5.30a) 



xgU^x + (^g + ,9^ Li = xgNi,^ + (^^g' - g^ Ni (5.306) 



xgXi^, 



-gLi 



-J-9]X, 



where 



Xi= Xi + Li, 



(5.30c) 



(5.31) 



which can be compared with the analogous Chan- 
drasekhar system ( [37], Ch.4, eqs. (46-47), (50)). Our 
equations (5.30b)-(5.30c) coincide with Chandrasekhar 
equations (46)-(47), while our equation (5.30a) coincides 
with the Chandrasekhar equation (50) if and only if 
p' = which is equivalent to {x'^g"/2 — g + 1) — 0. 

The basic system (5.30) can be directly applied to 
study extreme' black hole case. In the next subsection 
we investigate first the case of a simple horizon to make 
clear peculiarity of the case of the double horizon. 



B. Extreme black hole case 

Behavior near the simple horizon 

In the neighborhood of a simple horizon x+ we have 
g{x) = g'ix+){x - x+) + yix+){x - x+f + .... To 
study behavior in the limit a; ^ a;+ + we introduce the 
variable z = a: — a:+ . In a small neighborhood of ^; = 
limiting system for (5.30) reads 

x+{g'{x+))^z^Ni,, = {n + l)g'{x+)zNi 



+g'{x+)z 



x+ 



g'{x+)-{n+l) 



z{Ni-Li),, + -{Ni-Li)=Q 



Li+xl[-{g'{x+)fWijXi 
(5.32) 

(5.33) 
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-zLi + -^x+Xi 



(5.34) 



One immediately sees from (5.33) that the restricted 
near z = Q solutions should satisfy Ni{z) = Li{z). Then 
(5.32) and (5.34) form the system of two first-order equa- 
tions for functions Ni , Li 

{g'{x+)fz^Ni,, = ^{g'{x+fzNi+x+(\{g'{x+)f+af\xi, 



As follows from (5.42), for a restricted solution it should 
be Ni = Li. Then equations (5.41) and (5.43) form a 
system of two first-order equations for Ni,Li: 



x±-^Xi, 
T 



x+zXi^^ 



-zNi + ^x+Xi 



(5.35) 
(5.36) 



This system reduces to one second-order equation for Xi 



Z Xi^zz — zXi^^ 



1 + 



Xi=0 (5.37) 



This is the Euler equation, and solutions of (5.32)-(5.34) 
restricted near zero, have the form 



Xi{z) = 



Bii cos 



Ni=Li = -x+ 



In z +B21 sin 



7'(^+) 



In z 



z, 



(^+) 



B21 I cos 



(5.38) 



In; 



+ ( ^Ba; - I sin 

2 g'[x+) 



In; 



(5.39) 



where Bu, B21 are arbitrary constants. As a result in the 
small neighborhood of a simple horizon restricted solu- 
tions exist for all real values of ai . 

Behavior near the double horizon 

The double horizon x± corresponds to the case a = acr 
in (5.28). For the case of the density profile (1.4) 



1.456 



(5.40) 



In the small neighborhood of the point x = x±, the 
metric function is written as 

g{x) = 7(x - x±f + ...,7 = ]^g"{x±) 

In the variable z = x — x±, in the small neighborhood of 
z = Q the limiting system for (5.30) reads 

x±l^z^Ni^z = (n + l)^z\Ni - Li) + xlafXi, (5.41) 
ziNi - Li)^z + {Ni - Li) = 0, (5.42) 
x±zXi^z = -zLi + x±Xi (5.43) 



zXi^z = Ni+Xi 

x± 



(5.44) 



(5.45) 



This system reduces to one second-order equation for A^; 



z^Ni^zz + ^z^'Ni^z + %Ni = 



(5.46) 



which differs essentially from the analogous equation 
(5.37) for a simple horizon case. General solution to 
(5.46) is [44] 



Ni{z) 



CiiJi 



1 

7 z 



C21Y1 



1 

7 z 



(5.47) 



where Cu.C^i are arbitrary constants, Ji,Fi are Bessel 
functions. 

With taking into account asymptotic behavior of 
Bessel functions for big values of argument, we find the 
behavior of function Ni{z) for z ^ 



Ni{.z) = — 
z^ 



(ail 37r\ . / crj 1 37r\ 

Cii cos T + <-^2i sm - 

\^ z 4/ V7-^ ^ J 

(5.48) 

We see that solutions to (5.46) are unbounded as 2; ^ 
for all real values of the parameter cr;. 
From (5.44) we get 



Xi{z) = -- 



72 



CiiJo ( — - 
7 z 



{x±)ai 

which gives in the limit z 



7 z 



(5.49) 



Xi{z) 



{x±) \7r 



Cii cos - 

V 7 2: 4 





1 




(t 


z 





(5.50) 



Analysis of our basic system (5.30) in small neighbor- 
hood of double horizon x = x± shows that for all real 
values of the parameter ai, there exist unbounded solu- 
tions as X x±. Therefore the method of linear pertur- 
bations, as well as the assumption (5.8), are not suitable 
in this case, but the behavior of perturbations suggest 
instability of the extreme configuration. It should be in- 
vestigated separately, and we are currently working on 
this [45]. 
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C. The reduction of the system to a one-dimensional 
wave equation 

Now our goal is to reduce the system (5.30) to a single 
second-order equation. We introduce the new functions 
Zii, Z21, Z31 using the linear transformations 



N, 



g^, (5.51a) 



1 / 



x\2 

It is easily to prove that 
^x'g"' + xg"-g' 



g"' + xg"-g' ]z2i, 



(5.55c) 



— {x^p'Y = 3ax''p'^ (5.56) 



so that in the case when the density profile satisfies the 
condition 



X, 



g + TT 22; + Z31 



(n + 1) g X ., 

--- + 7^9 \Z21-Z31 

X X 2 



(5.516) 



(5.51c) 



The inverse transformation to (5.51) reads 



ziiix) 



V9b 



b{x)Ni{x)+{{b{x)-n-l) 



xg_ 
2 9 



xb' 



-'■<-)+<''i&)^'+(('-"-47-''-+"''^)* 

(5.52a) 



Z2i{x) = [Li + Xi] 



(5.526) 



(xV)' = 0, 



(5.57a) 



the system (5.55) splits on the system of two equations 
(5.55a), (5.55c) for z-^uz^i, and the equation (5.55b). 
Condition (5.57a) is in turn equivalent to 



(5.576) 



In the particular case [x^p') = const = this is 
the necessary and sufficient condition for coinciding of 
OTir system (5.30) with the Chandrasekhar system ( [37], 
cqs.(46)-(47), (50) Ch.4). 

Differentiating (5.55c), we come to the system which 
includes one first-order equation, (5.55b), and one 
second-order equation 

Z3i,xx + ^ (^^ + '^3!,x + qi {x)z3i = n {x)z2i , (5.58) 
where 



Z3l{x) 

where 



(5.52c) 



b{x) =n+l + -g\x) - g{x) = « + ^ ( M{x) - x^p 



The sum of (5.51b) and (5.51c) gives 
Xi + Li = 



n + l 1 , 1 

-ir^2'--x' 



(5.53) 



g^l^Z2i (5.54) 



As a result we get the following system 



x' fg" (gT 2I 



b{x) V 2 4g 



Z31, (5.55a) 



2 1 2(n + l) Ig" , Ifg' 



qi{x)=ai^- 



2 g 4\gJ x g 



{xg"-g') 



b{x) 



{xg"-g') g' ^ 1 
26(x) g X 



xg 



iax 



2b{x) ' b{xf^ 
(5.59) 



ri{x) = -3apY 



jn+l) _ 3xg^ 
~9 ~^~9 



+ ^(xV^)'. 



(5.60) 



It is easily to sec that in the case when the condition 
(5.57) is satisfied, two equations (5.55b) and (5.58) split. 
Introducing the new function u)3i{x) by 



Z3l{x) = —0J3l{x) 

xg 



(5.61) 



Z2l,x = — 



b{x) 



(5.556) 



Z3L: 



b{x)g- 



-zii 



2 ^ (xg" - g') 

X 



2b{x) 



Z31 



we reduce the equation (5.58) to the form which does not 
contain the first derivative: 



1^31,: 



'^f^-Vnix) 



W3i = xgriZ2i, 



(5.62) 
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where the potential Vu (x) is given by 

1(1 + 1) 1 , 3g" Ig' , ixg"-g') 

Vll[X) = 2 To 1 U~\ 

g 2 g X g b{x) 



{xg"-g') 
2b{x) 



5 xj A\gJ 2b{x) xbix) [2'" ^ ^ 

(5.63) 

With taking into account (5.55b) and (5.61), equation 
(5.62) can be rewritten as integro-differential equation of 
the form 



g-\x)af -Vu{x) 



(jJaiix) 



-xg{x)ri{x) j 



xg{x)b{x 



■dx 



(5.64) 



In this form {g~'^{x) scales the spectral parameter a^) the 
equation (5.64) corresponds to the generalized spectral 
problem with the non-local potential 

-W3;,xx + Vu{x)u!si{x) - Tiu}3i{x) = af\uj3i{x), (5.65) 

9 

where 

X 

u{z)dz 



Tiu(x) = xg{x)ri{x) 



zb{z)g{z) 



(5.66) 



is the integral Vol'terra operator. The lower limit is = 
x+ for a black hole case. 

The condition (5.57) (which leads to r; = 0) is nec- 
essary and sufficient condition to reduce (5.64) to the 
Schrodingcr equation with the local potential. 

Introducing "the tortoise coordinate" a;* (a;) = 
J dx/g{x) and the function w;(.t») by 



W3i{X:,) = x^/g{x)z3i{x*) (5.67) 
we reduce the system (5.58), (5.55b) to the form 

W3i,x^x, + - Wi{x)]w3i{x^) = xgi{x)ri{x)z2i{x^) 

(5.68) 

(5.69) 



,1/2 



where 

Wi{x) 



1(1 + , 



xb{x) 



1 , g{j:g"-g') ^ a:g"-g' 

x^ b \ 2b 



xgg 
2b 



xb\2 



x'g'" + xg"-g' 



(5.70) 



In the limit x the integral term in (5.66) tends to 

zero, on essence due to Z21 0. Indeed, when x ^ x+, 
we get Z21 ^ ^yx — Xj^ with using (5.52b) and taking into 
account asymptotic behavior of (L; + Xi) ^ {x — x^) 
which follows from (5.38). 

The potential (5.70) vanishes as a;* — > -|-oo as x~'^, 
while for —00, it vanishes exponentially. Therefore 



solutions to (5.68) have asymptotic e 



as x^ 



±00, 



so that we have to look for solutions satisfying boundary 
conditions 



Wl 



CIS t^s^i 



00 



Wl ^ r/^'^e* 



as a;* —00 



(5.71) 



In the particular case of validity of (5.57), ri{x) = 0, 
p'j_ = 0, the system (5.68)-(5.69) splits and we get the 
Schrodinger equation 

-W3i,x,x, + Woi{x)w3i{x^,) = afw3i{x^.) 
with the potential 



Woi{x)=g 



x^ 



+ 



gixg" - 9') ( xg" - g' 



b{x) 



2b{x) 



9_ 

9 



X 

2 

X 



which for the Schwarzschild geometry coincides with the 
potential in the Zerilli equation ( [37], Ch.4, eq.(63)). 

We have reduced the basic system of three first-order 
linear equations (5.29) to a single second-order equation 
for the particular combination of these function, W3;(a;»). 
In our case this is the Schrodinger equation (5.72) with 
non-local potential. Its non-local part vanishes when the 
condition (5.57) is satisfied. As in Schwarzschild case, 
this empirically found reducibility (resulted from a se- 
quence of mysterious cancellations), follows in fact from 
the existence of some particular solution to the system 
which we have to reduce [37]. 

In our case the condition (5.57) is the necessary and 
sufficient condition for vanishing non-local part in (5.72) 
and thus for reducing our system (5.29) to the standard 
Schrodinger equation. Actually, the condition (5.57) is 
also the necessary and sufficient condition for existence 
of particular solution which guarantees such a reduction 
(it is shown in detail in our paper [45] devoted to the case 
satisfying (5.57)). If the condition (5.57) is satisfied, then 
the particular solution reads 



N! = V9 



2g 



29 



-at-^-b'{x) 



L\ = -^gb'{x); Xf = ^ 



b{x) 



b'ix) 
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The linear transformation (5.51) is needed to reduce 
the system of 3 hnear equations to one equation, if the 
particular solution to (5.29) exists. Remarkable luck in 
our case is that this transformation works not only in 
the case when (5.57) is satisfied (i.e. obtained second 
order equation is the standard Schrodinger equation with 
the local potential), but also in general case. This fact 
allowed us to reduce our system (5.29) to the Schrodinger 
equation (5.72) with the non-local potential. 

We see that the problem of stability of the spherically 
symmetric solution (2.7) to polar perturbations, reduces 
to investigation of the spectral problem (5.68)-(5.69) with 
the potential of the form (5.70) and with the boundary 
conditions (5.71). 

Equation (5.68), with taking into account (5.69), can 
be written in the form corresponding to the spectral prob- 
lem with a non-local potential 



(5.72) 



where 



fiu{x^) = x{x^)g^{x{x: 



■■))ri{x{x.)) J 



g'^{x{z^))u{z^)dz^ 
x{zt,)b{x{Z:t)) 

(5.73) 

is the integral Vol'tcrra operator. 

Ultimately the case reduces to investigation of the 
eigenvalue problem with integro- differential operator 
(5.72), which is the Schrodinger equation with the non- 
local potential. 

The spectrum of eigenvalues contain all the values of 
the parameter af, at which solutions exist which satisfy 
the imposed boundary conditions. 

If such solutions exist only for real values of the time 
parameter ai and if, in addition, they form the complete 
basic set of functions, then any smooth initial perturba- 
tion on the finite interval of variable x» (with compact 
measure) can be expanded on this basic set, and since 
dependence of each particular mode on time is given by 
exp(i(T;f), it testifies for stability of geometry. 

Indeed, a considered static configuration is stable if 
there are no intcgrable modes with negative af . Appear- 
ance of negative eigenvalues af would lead to existence 
of exponentially growing modes of perturbations. 

In the next subsection we study in detail the integro- 
differential operator governing polar perturbations. 



D. Schrodinger equation with non-local potential 

A system governed by Schrodinger equation with a 
non-local potential, obeys the following theorems: 

i) If in the standard one-dimensional Schrodinger equa- 
tion the potential is nonnegative, then negative eigenval- 
ues are absent (see, e.g., [42]). 



ii) The Weyl theorem [46] for self-conjugate operators: 
The essential spectrum conserves under relatively com- 
pact perturbations [47]. 

The essential spectrum is defined as follows: If we re- 
move from the spectrum of self-conjugate operator all 
isolated points which are eigenvalues of finite multiplic- 
ity, then the remaining of the spectrum is the essential 
spectrum. 

Essential spectrum of non-perturbed (local) poten- 
tial is continuous and represented by positive semi-axis 
[0, oo), isolated points are absent in case when negative 
values are excluded by non-negativity of the potential. 

It follows that the essential spectrum of the problem 
with the non-local potential (5.72) is the same as the 
essential spectrum of the non-perturbed (local) poten- 
tial. Essential spectrum of local potential is this total 
spectrum, since isolated points are absent, i.e., essential 
spectrum of the perturbed problem coincides with the 
total spectrum of the non-perturbed problem. 

Non-local part of a potential represents the perturba- 
tion of the local potential. To not spoil an essential spec- 
trum, this perturbation should be relatively compact. 

So, our task now is to prove that non-local part repre- 
sents a compact perturbation and to deduce criterion of 
non-negativity of a local potential. 

In our case non-local part (perturbation of a local po- 
tential) is given by the integral Vol'terra operator (5.73). 
Such an operator is totally continuous, if it has smooth 
square integrable kernel. 

Square integrability requires 



/oo PC 
-oo J —( 



K'^{Xt:, Z:i)dX^dZt: < OO 



(5.74) 



The kernel of our Vol'terra operator (5.73) is 



K{x^,,Zt,) = x{x*)g^ {x{x*))ri{x{x*))- 



9^{x{z*)) 



x{zt:)b{x{zt,)) ' 

(5.75) 

Its smoothness is evident. The sufficient condition for 
square integrability of the kernel K{x,y) is the condition 
on behavior of K'^ at infinity: 



K\x,y) < 



1 



"l + i5i ^,1+62 



where 5i, S2 are arbitrarily small. 

For K{x^., z,) given by (5.75), for —00, K"^ van- 

ishes as (7(x(z,)). When —00, then K"^ vanishes 

as g^{x{xt,)). The metric function g{x{Xi,)) near hori- 
zon behaves as g{xt) ~ g'{x^)e^ (x,)x,^ gi jg positive, so 
that metric as a function of a;* vanishes exponentially at 
approaching the horizon. 

When z* ^ 00, then K"^ vanishes as a;~^(z»). Prom 
definition z» we see the main contribution at infinity 
is 2;* ^ X. When a;* — > 00, K"^ vanishes as a:^(p^)^. 
The tangential pressure p_L for de Sitter-Schwarzschild 



14 



geometry vanishes at infinity quicker than , because 
it is related with density by the equation of state p± = 
—p — xp'/2, and density vanishes quicker than to 
guarantee the finitcness of a mass. Hence in this Umit 
vanishes quicker than x~^. 

So, for a BH case the kernel is square integrable. 

As a result the totally continuous operator (5.73) gives 
a relatively compact perturbation to the local potential 
in the integro- differential equation (5.72). 

Criterion of non-negativity of local potential 

Introducing the function p{x) = xg"{x) — g'{x), we 
write the potential (5.72) in the form 



Wi{x) = g 



where 



1 



2h{x) 



ig') 



2(n+ 1) 



1 

bx 



(5.76) 



Il{x)=x'[^g'g" + gg"' 



-{n + l-g)p{x) (5.77) 



In (5.76) we should investigate the term J;. The rest is 
positive, since b{x) > n in a BH case. 

Expressing g{x), its derivatives, p{x) and b{x) in terms 
of mass function M{x), density p{x) and its derivatives, 
we transform Ii{x) to the form 



Ii{x) 



-4a- 



M2 



+ -ax^pp' + 3(n — l)x^p' 



> 



> ^(n-2) >0. 



(5.80) 



As a result, we find the sufficient condition (5.78) for 
non-negativity of the potential (5.70) in all range of ar- 
gument for which g{x) > 0. 

For the density profile (1.3) this condition is satisfied. 

We can conclude that the essential spectrum of the 
integro-diffcrcntial operator (5.72) is the same as the es- 
sential spectrum of its local potential. Now the key point 
is to find the condition on a perturbation of a local po- 
tential which guarantees the absence of isolated points in 
the total spectrum (negative values of a'f ) of the integro- 
differential operator (5.72). 



Non-local contribution 

Multiplying (5.68) by wt,i* and integrating by parts 
with taking into account asymptotic behavior of (5.60) 
at infinity, we obtain the following relation 



+00 

(^f y \wu{x*)\'^dx*-\-w:ii^x,wli\ 



+ 00 

/" 



[I'U'si.xJ^ + Wi{x)\wzi{x^)\^ + i^i,x9{x)\z2i{x^)\^]dx^. 



(5.81) 



where 



+ -axp'M - Snxp - 3gx(x^p" -\- 2p) + ^'■"'^^W 
2 x-^ 

(5.77a) 

For a BH case g{x) > Q while the weak energy condition 
gives p' < 0. Then the sufficient condition for W; > is 
the condition on the equation of state 

x^p"{x) + 2p{x) > 0, (5.78a) 

which constraints the growth of the derivative oip_\_-\- p 

x{p± +p)' <P+ {p± + p) (5.786) 

This condition is actually satisfied also for the case with- 
out horizons (then g{x) > Q for all a;). 

When (5.78) is satisfied, then proof of non-negativity of 
(5.76) reduces to proof of non-negativity of the function 



cl>{x) = 



2(n-hl) ia{n + 2)M 



x'^b{x) 

It is bounded from below as follows 



(5.79) 



cl>{x) = -3 



n+l 



3aM{n + 2) 
2xb{x) 



i^ii^) = -^9^{x)b{x)ri{x) 



(5.82) 



The Wronskian Wzi,x,wli\^'^ of two independent solu- 
tions and W3; is constant (see [37], Par. 27 Ch.4). 

The contribution to the spectrum from the non-local 
part of the potential is given by 



-t-CX) 

N= j i}i,x{x)\z2i{x)\'^dx 

x+ 



(5.83) 



If the condition of non-negativity of a local potential 
Wi{x) is satisfied, then the requirement 

+00 

N= j i^iAx)\z2i{x)\''dx>Q (5.84) 



gives the sufficient condition for the absence of negative 
eigenvalues of the considered spectral problem. 



2 

> — 



, 2,{n + 2)a.M' 

n+1--^ ^ 

Z n X 



, 3(n + 2),_, , 
n+l--^(l-,) 



*We denote the complex conjugate by * for convenience of 
comparison with the classical results presented in [37]. 
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Fortunately, non-local contribution given by (5.83), 
does not grow with the mode number n, since \z2i\'^ is 
constrained from above by the function proportional to 
n~^. This constraint is valid for any density profile and 
follows from (5.69) with taking into account that in a BH 
case h{x) > n. 

In the case when the metric function g{x) satisfies the 
condition (5.57), the sufficient condition (5.84) is triv- 
ially satisfied (iV = 0), and negative eigenvalues do not 
appear in the spectrum. As a result spherically symmet- 
ric metrics satisfying (5.57) and (5.78) are stable to polar 
perturbations. 

VI. THE CASE OF DENSITY PROFILE (1.3) 

In the case of the density profile p{x) = e~^^ the metric 
fimction in (5.28) reads 

5(x) ^e^'^^") = 1- -(1-e-"') (6.1) 

X 

Potential (4.14) governing the axial perturbations is 
depicted in Figs. 5-6 for the density profile (1.3) and two 
values of the characteristic parameter a (denoted in fig- 
ures as a). 
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Potentials are smooth, short-range and positive for all 
values of the characteristic parameter a. Therefore all 
types of vacuum configurations with de Sitter center in- 
cluding a vacuum nonsingular black hole, are stable to 
axial perturbations. 

The local potential governing polar perturbations 
given by (5.70) is shown in figs. 7-8 for the density profile 
(1.3) and two values of parameter a (denoted in figures 
as a). 
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FIG. 7. Polar potential (5.70) for m 2.Smcrit- 
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FIG. 5. Axial potential (4.14) for m ~ 2.8mc 



FIG. 8. Polar potential (5.70) for m ~ 11.2mc 
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FIG. 6. Axial potential (4.14) for m ~ 11.2mc 



Both axial and polar potentials for bigger values of a 
become similar to those for the Schwarzschild case [37] . 

Polar local potentials are smooth short-range poten- 
tials, so that integrals of them are finite over all the re- 
gion of variable x. 

The potential (5.70) for the density profile (1.3) satis- 
fies the criterion of non-negativity (5.78), but the condi- 
tion (5.57) is not satisfied, so that appearance of negative 
eigenvalues af is in principle possible. 



The question of existence of isolated points with neg- 
ative values af in the total spectrum of the integro- 
differential operator(5.72) requires the complicated nu- 
merical analysis which is in progress. Preliminary results 
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suggest that non-local contribution (5.83) does not lead 
to negative values for the masses m > mcr2- 

This result looks natural. The second critical mass 
value TOcr2 is distinguished for the unperturbed geometry. 
The value rricr^ marks the point in the temperature-mass 
diagram at which specific heat is broken and changes 
its sign, so that a second order phase transition starts 
when in the course of Hawking evaporation the mass ap- 
proaches the value mcr2- For the density profile (1.3) it 
is given by 

mcr2 ^ 0.38mpi^yppi/pQ (6.2) 

The extreme state of non-singular black hole (m = rricr), 
can be unstable since some perturbations modes grow 
unlimited at the double horizon for any density profile. 
If the considered configuration would develop instability 
before achieving the extreme state, the most appropriate 
range gets beyond mcr2 where a phase transition starts. 

The critical value mcr2 corresponds to the maximum 
at the temperature-mass curve (see fig. 3). It is calculated 
from the condition dT/dm = 0. In the units normalized 
to de Sitter radius ro (which is the characteristic scale 
related to de Sitter vacuum trapped in the origin) 

y+ = ^; (6.3) 
the temperature on a BH event horizon is given by [17] 
T=±-±(l-y±] (6.4) 

y+ y+\ s J 

The density profile and metric in these units read 

p{y) = e-y'/^; g{y) = l-i(^l-e-y'/'^ (6.5) 

Prom dT/ds = and g{y+) = 0, we get the critical value 
S2 and the value y+ corresponding to m = rricri- 

S2 ==i 2.226; y+ =~ 2.166 (6.6) 

For comparison, the critical values for the extreme case 
of the double horizon r± are [17] 

Scr ^ 1.7576; y± ~ 1.4957 (6.7) 

VII. DISCUSSION 

We present the conditions specifying two types of con- 
figurations with de Sitter center, including black holes 
with and without changes of topology. 

We found that any configuration described by spheri- 
cally symmetric geometry with de Sitter center is stable 
to axial perturbations. 

The problem of stability to polar perturbation reduces 
to a one-dimensional Schrodinger equation with a non- 
local potential given by the Vol'terra integral operator 



with square integrable smooth kernel representing a com- 
pact perturbation to the local potential. 

We derived the criterion of non-negativity of the local 
potential which defines the essential spectrum of integro- 
differential operator governing polar perturbations in 
general case. 

We derived the criterion of vanishing of non-local part 
of the potential which distinguishes the class of geome- 
tries for which the problem of stability reduces to a stan- 
dard one-dimensional Schrodinger equation. This class 
of metrics has been studied in [45] . 

For the case when perturbations are described by the 
Schrodinger equation with the non-local potential, we 
found the sufficient condition for the absence of the neg- 
ative eigenvalues in the spectrum which guarantees the 
stability of investigated geometry. 

For an extreme black hole, the method of small per- 
turbations is not applicable due to existence of unlimited 
perturbation modes at approaching the double horizon. 

Asymptotic behavior of the basic system near double 
horizon suggests instability of the extreme configuration. 
The behavior in this regime is very special, unrestricted 
solutions for perturbations exist for positive values of the 
spectral parameter af. The limiting equation for pertur- 
bations near double horizon is essentially different from 
that for one-horizon case in which unrestricted solutions 
do not appear for positive values af and which cannot 
be smoothly continued to two-horizon case. The ques- 
tion arises what is the place of the metric with m = rricr, 
at the set of metrics whose stability we investigate as 
one-parametric set of solutions to Einstein unperturbed 
equations with a given density profile. 

The critical value of the mass parameter rricr is cal- 
culated from two transcendental equations: g{r±) = 
0; 9'{f±) — 0. This is the unique point for considered 
metric function g{r) because it is the minimum of g{r) 
and g{r) has the only one minimum [20]. Therefore the 
transcendental system for {r±,'mcr) has the unique solu- 
tion for each particular one-parametric set with a given 
density profile. 

The metric with double horizon represents an isolated 
singular point at the set of metrics g{r) for each given 
density profile. It resembles a fixed point attractor be- 
havior which is currently the key point of our efforts [45] . 
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